Quantum Simulation of the Majorana Equation and Unphysical Operations 
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We design a quantum simulator for the Majorana equation, a non-Hamiltonian relativistic wave 
equation that might describe neutrinos and other exotic particles beyond the standard model. The 
simulation demands the implementation of charge conjugation, an unphysical operation that opens 
a new front in quantum simulations, including other discrete symmetries as complex conjugation 
and time reversal. Furthermore, we describe how to implement this general method in trapped ions. 



A quantum simulator is a device engineered to repro- 
duce the properties of an ideal quantum model. This 
still-emerging topical area has generated a remarkable 
exchange of scientific knowledge between apparently un- 
connected subfields of physics. In terms of applications, 
it allows for the study of quantum systems that can- 
not be efficiently simulated on classical computers pQ. 
While a quantum computer would also implement a uni- 
versal quantum simulator [2], only particular systems 
have been simulated up to now using dedicated quan- 
tum simulators [3]. Still, there is a wealth of successful 
cases, such as spin models [H[5], quantum chemistry [5] 
and quantum phase transitions [7]. The quantum simu- 
lation of fermionic systems [8] and relativistic quantum 
physics have also attracted recent attention, reproduc- 
ing dynamics and effects currently out of experimental 
reach. Examples include black holes in Bose-Einstein 
condensates [5] , quantum field theories [TU1 |TT] and re- 
cent quantum simulations of relativistic quantum effects 
as Zitterbewegung, Klein paradox and interacting rela- 
tivistic particles [T^HTo] in trapped ions. 

In this paper, we show how the Majorana equation [T7] 
can be simulated in an analog quantum simulator, hav- 
ing as a key requirement the implementation of complex 
conjugation of the wavefunction. In this manner, we are 
able to propose this and other unphysical operations such 
as charge conjugation and time reversal, constituting a 
novel toolbox of accessible quantum operations in the 
general frame of quantum simulations. While quantum 
simulators may soon realize calculations impossible for 
classical computers, we show here the possibility of im- 
plementing quantum dynamics that are impossible for 
our quantum world. 

The Majorana equation is a relativistic wave equation 
for fermions where the mass term contains the charge 
conjugate of the spinor, ip c , 

ihtfip = mcip c . (1) 

Here, $ = j^d^ and 7 M are the Dirac matrices jTHj, while 



the non-Hamiltonian character stems from the simulta- 
neous presence of ip and yb c . The significance of the Ma- 
jorana equation rests on the fact that it can be derived 
from first principles in a similar fashion as the Dirac equa- 
tion [171 119] . Both wave equations are Lorentz invariant 
but the former preserves helicity and does not admit sta- 
tionary solutions. The Majorana equation is considered 
a possible model [20] for describing exotic particles in 
supersymmetric theories -photinos and gluinos-, or in 
grand unified theories, as is the case of neutrinos. In- 
deed, the discussion of whether neutrinos are Dirac or 
Majorana particles still remains open [2T]. Nevertheless, 
despite the similar naming, this work is neither related 
to the Majorana fermions (modes) in many-body sys- 
tems [221 H3] j uor to the Majorana fermions (spinors) in 
the Dirac equation [2"D"ll2"3] . 

In order to simulate the Majorana equation, we have to 
solve a fundamental problem: the physical implementa- 
tion of antilinear and antiunitary operations in a quan- 
tum simulator. Here, we introduce a mapping [25] by 
which complex conjugation, an unphysical operation, be- 
comes a unitary operation acting on an enlarged Hilbert 
space. The mapping works in arbitrary dimensions and 
can be immediately applied on different experimental se- 
tups. We show how to simulate the Majorana equation 
in 1+1 dimensions and other unphysical operations us- 
ing two trapped ions. We also give a recipe for measuring 
observables and a roadmap towards more general scenar- 
ios. In this sense, this work provides a novel toolbox for 
quantum simulations. 

There are three discrete symmetries [26] which are cen- 
tral to quantum mechanics and our understanding of 
particles, fields and their interactions: parity, V, time 
reversal, T, and charge conjugation, C. None of these 
operations can be carried out in the real world: V in- 
volves a global change of physical space, while C and 
T are antiunitaries. However, there is no apparent fun- 
damental restriction for implementing them in physical 
systems that simulate quantum mechanics. We will fo- 
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FIG. 1. Diagram showing the different steps involved in the quantum simulation of unphysical operations in 1 + 1 dimensions. 



cus on the study of antiunitary operations, which can be 
decomposed into a product of a unitary, Uc or U-y, and 
complex conjugation, Kip — %p* . We consider the map- 
ping of the quantum states of an n-dimensional complex 
Hilbert space, C w , onto a real Hilbert space [25], M.2 n , 



<p) 



(2) 



This mapping can be implemented by means of an auxil- 
iary two- level system, such that Bbn G ^2 <£> Hn . In this 
manner, the complex conjugation of the simulated state 
becomes a local unitary Vk: acting solely on the ancillary 
space, Kip — i]j* — » = (a z ®l) i &, and thus physically 
implementable for a wavefunction of arbitrary dimen- 
sions. Furthermore, unitaries and observables can also be 
mapped onto the real space, 0—^0 = 1® O r — io v ® Oi, 



where O r = \{0 + KOK) and O, 



r(0 - KOK), 



preserving unitarity and Hermiticity. The proposed sim- 
ulator also accommodates the antiunitary operations C = 
IAqK and T = Uj-K. To this end, we have to choose a 
particular representation that fixes the unitaries Uc and 
U-y, as will be shown below. 

We possess now the basic tools to simulate the Majo- 
rana equation (JlJ. The expression for the charge con- 
jugate spinor is given by ip c = W7 Kip, with W a uni- 
tary matrix satisfying W~ 1, ~f tJ "W = — (7 P ) T ■ We illus- 
trate now the proposed quantum simulation with the case 
of 1+1 dimensions. Here, a suitable representation of 
charge conjugation is ip c = ia y a z ip* , that is W = icr y , 
and the Majorana equation reads 



ihdfip = ccr x p x ip — imc 2 <jyip* 



(3) 



where p x — —ihd x is the momentum operator. Note that 
Eq. Q is not Hamiltonian, (ihdtip 7^ Hip). This is due 
to the presence of a complex conjugate operation in the 
right-hand side of Eq. ([3]), which is not a linear Hermitian 
operator. Surprisingly, through our mapping 



(4) 



the Majorana equation for a complex spinor becomes a 
3+1 Dirac equation with dimensional reduction, p y ,p z — 
0, and a four-component real bispinor 



ihdt^ = [c(l ® cr x )p x 



mc 2 <7 x 



(5) 



Here, the dynamics preserves the reality of the bispinor 
^> and cannot be reduced to a single 1+1 Dirac particle. 
In general, the complex-to-real map in arbitrary dimen- 
sions transforms a Majorana equation into a higher di- 
mensional Dirac equation |27j . Since Eq. § is a Hamilto- 
nian equation, it can be simulated in a quantum system. 

The mapping of wavefunctions into larger spinors also 
allows us to explore exotic symmetries and unphysical op- 
erations, otherwise impossible in nature. From Eqs. pj), 
Q, and |5|, for the 1+1 dimensional case, we can deduce 
that charge conjugation is implemented in the enlarged 
space via the unitary operation Vq 

<p c = Cip= U c Kip -> V c 9 = -(<r g ® <J X )K- (6) 

We can do something similar with time reversal, defined 
as the change t — > (— t). In this case, we expect [T§] 
ihd T ip' [t) — HiP'(t), where the time variable r = —t 
and the modified spinor ip'(r) = Tip(t). In order to pre- 
serve scalar products and distances, the time reversal 
operator must be an anti-unitary operator and thus de- 
composable as the product T = Uq-K. In 1 + 1 dimen- 
sions, imposing that the Hamiltonian be invariant under 
time reversal, T~ X HT , implies that the unitary satisfies 
U^ 1 {i<j x d x )lAj- — —ia x d x , with a possible choice being 
Uj- = a z . In other words, in the enlarged simulation space 



Tip = U T Klp -+ = ( a z ® crj*. 



(7) 



See Fig. [T] for a scheme of the simulated symmetries. As 
mentioned before, quantum simulations of unphysical op- 
erations can be straightforwardly extended to higher di- 
mensions. In this sense, Eqs. ^ and ^ will be valid for 
wave functions ip of dimension d as long as we consider 
the complex conjugation of an arbitrary wavefunction as 
= (a z ® Id)*. 
The proposed protocol for implementing unphysical 
operations onto a physical setup allows us to deal with 
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situations that are, otherwise, intractable wih conven- 
tional quantum simulations. To exemplify the value of 
this novel building block in the quantum simulation tool- 
box, we consider the case of an advanced experimental 
quantum simulation, impossible to reproduce with classi- 
cal computers. We assume that, after a certain evolution 
time, it is crucial to realize an unphysical operation such 
as charge conjugate or time reversal, before continuing 
the unitary (physical) evolution. With existing tools in 
quantum simulations, we would need to stop the dynam- 
ics, implement a full quantum tomography of the current 
quantum state associated to a huge Hilbert space, apply 
the unphysical operation in a classical computer, encode 
back the modified quantum state into the experimental 
setup, and then to go ahead with the quantum simula- 
tion. Clearly, this task would be impossible with classical 
resources and would become possible with a suitable im- 
plementation of our proposed ideas. 

In a recent experiment, the dynamics of a free Dirac 
particle was simulated with a single trapped ion |13j . 
Here, Eq. ^ has a more complex structure requiring 
a different setup. Moreover, the encoded Majorana dy- 
namics requires a systematic decoding via a suitable re- 
verse mapping of observables. We can simulate Eq. ^ 
in two trapped ions, with lasers coupling their internal 
states and motional degrees of freedom. The kinetic 
part, cp x (t ® cr x ), is created with a laser tuned to the 
blue and red motional sidebands of an electronic transi- 
tion [T2l [14] , focussed on ion 2. The second term, a x ®a y , 
is derived from detuned red and blue sideband excitations 
acting on each ion. The Hamiltonian describing this sit- 
uation reads 



H = h-^-o\ + h-^-a^ + hva) a + hu r b^b 

+ nn r( e i(9*i-'"i*+*i) + e <(9*i-<4*+*i))< 7 + + h 

+ Ml Ue^-^t+M + e i(qZ2-u,' 2 t+j,' 2 )j a + + jj 
+ Ml [( e *(8*-^+(t) + e i(.qZ2-u't+<t>')j a + + H c 



Here Zip = Z ± § are the ion positions, measured from 
the center of mass, Z, and relative coordinate, z. The 
phases of the lasers fa for i = 1,2, (fa fa), are controlled 
to perform the interaction term (kinetic term). The fre- 
quencies of the center of mass and stretch mode are given 
by v and v r = \[?>v, while a) , a, b\ and b, are the corre- 
sponding creation and annihilation operators. Finally, £1 
and are the laser Rabi frequencies in the rotating-wave 
approximation. With the adequate choice of parameters, 
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(8) 



the Hamitonian in the interaction picture reads 

H = Hr] r Q(a x <g 1 - 1 ® <r y )(tfe iSt + be~ iSt ), 
+hrjQ(TL <g <j x )i(a) — a) 



(9) 



where r\ = ?7 r 3 1 / 4 = \J h/Am'v « 1 is the Lamb- 
Dicke parameter and m' the ion mass. In the limit of 
large detuning, we have S ^> r] r Cl^ (frffr), T)Q\{al — a)\ 
and we recover Eq. ^ with the momentum operator 
p x — ih((v— a)/2A and the equivalences c = 2ryAf2 and 



mc 2 = 2hr) 2 .Q 2 /6 with A = ^h/Am'v. Introducing the 
ratio 7 = |mc 2 / (cp x )\, with 7 — gfajVj) 



- , we see it 



|(i( ra t-o)>|(77fi/5) 5 

is possible to tune the numerator and denominator inde- 
pendently so as to preserve the dispersive regime, while 
exploring simultaneously the range from 7 ~ (ultrarel- 
ativistic limit) to 7 — > 00 (nonrelativistic limit). 

A relevant feature of the Majorana equation in 3+1 di- 
mensions is the conservation of helicity. A reminiscent of 
the latter in 1 + 1 dimensions is the observable called here- 
after as pseudo-helicity £ = cr x p x . This quantity is con- 
served in the 1+1 Majorana dynamics of Eq. ^ but not 
in the 1+1 Dirac equation. We will use this observable 
to illustrate measurements on the Majorana wavefunc- 
tion. The mapping for operators can be simplified if we 
are only interested in expectation values. Reconstructing 
the complex spinor with the non-square matrix -0 = 
and M — (l it ), associated with Eqs. Q and (JHj), we 
have (0) 4 , = (ip\0\ip) = (^\M^OM\^) =: (O)*. There- 
fore, to obtain the pseudo-helicity £, we have to measure 



£ = M^a x p x M = (1 ® a x - a y ® a x ) ® p x 



(10) 



in the enlarged simulation space. In ion-trap experi- 
ments, we can use laser pulses to map information about 
the pseudo-helicity onto the internal states. The appli- 
cation of a state-dependent displacement operation on 
ion 2, U2 = cxp(— ik(t (g a y ) ®p x /2), generated by res- 
onant blue and red sidebands, followed by a measure- 
ment of 1 (g <r z , is equivalent to measuring the observ- 
able [13] A(k) = cos(kp x )('t (g a z ) + sm(kp x )(t (g a x ). 
Here, k is proportional to the probe time t pro b e . Note 
k=o ^ {{^-® a x)®Px)- Therefore, the first 



that ±(A(k)) 



term in ( 10 ) can be measured by applying a short probe 



d , 



pulse to the ions and measuring the initial slope of the 
observable A(k). To measure the second term in Eq. ( |To| , 
we apply the operation U\ — exp(—ik(a x <g 1) ®p x /2), 
and measure the spin correlation a z ®a x . We have, then, 
Vx)\ k=0 = 2((a y <g a x ) ®p x ). 
So far, we have presented a complete toolbox of un- 
physical operations, C, T, and JC. We can combine all 
these tools to study dynamical properties of the trans- 
formed wavefunctions. To exemplify the kind of experi- 
ments that become available, we have studied the scat- 
tering of wavepackets against a linearly growing poten- 
tial, V(x) = ax, with conventional numerical tools. It 
is known that repulsive potentials can be penetrated by 
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FIG. 2. Scattering of a fermion against a linearly growing potential (inset), (a) Ordinary Klein process (b) At an instant of 
time we apply the time reversal operator T causing the particle to retrace its own trajectory, (c) Similar to (b) but now we 
apply charge conjugation, converting the particle in its antiparticle. (d) Scattering of a Majorana particle, which propagates 
through the potential. Parameters are m = 0.5, c = 1 and V(x) — x, in dimensionless units. 



Dirac particles [18] , due to Klein tunneling [21 [15] . This 
is shown in Fig. [2^,, where a Dirac particle splits into a 
fraction of a particle, that bounces back, and a large an- 
tiparticle component that penetrates the barrier. This 
numerical experiment has been combined with the dis- 
crete symmetries and the Majorana equation. In Fig. [2]3 
we apply the time reversal operation after the particle 
has entered the barrier: all momenta are reversed and 
the wavepacket is refocused, tracing back exactly its tra- 
jectory. In Fig. [2]j we apply charge conjugation, chang- 
ing the sign of the charge and turning a repulsive elec- 
tric potential into an attractive one, which can be easily 
penetrated by the antiparticle. In Fig. [2]i, we show the 
scattering of a Majorana particle. While there are no 
plane wave solutions in the Majorana equation, we can 
still see a wavepacket penetrating the barrier, showing a 
counter-intuitive insensitivity to the presence of it. 

The previous implementation of discrete symmetries is 
valid both for Majorana and Dirac equations. Equally 
interesting is the possibility of combining both Dirac and 
Majorana mass terms in the same equation |20j . iH$?p = 
rriMCtpc + mDCip, which still requires only two ions. It 
also becomes feasible to have CP violating phases in the 
Dirac mass term, mrj exp(i#7 5 ). Furthermore, we could 
study the dynamics of coupled Majorana neutrinos with 
a term Mip c , where M is now a matrix and tp = ip(xi , X2) 
is the combination of two such particles, simulated with 
three ions and two vibrational modes. 

In summary, we have introduced a general method to 
implement quantum simulations of unphysical operations 
and non-Hamiltonian dynamics, such as the Majorana 
equation, in a Hamiltonian system. 
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